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ABSTRACT 

NGC 3201 is a globular cluster suspected to have an intrinsic spread in the iron content. We re-analysed a 
sample of 21 cluster stars observed with UVES-FLAMES at the Very Large Telescope and for which Simmerer 
et al. found a 0.4 dex wide [Fe/H] distribution with a metal-poor tail. We confirmed that when spectroscopic 
gravities are adopted, the derived [Fe/H] distribution spans ~0.4 dex. On the other hand, when photometric 
gravities are used, the metallicity distribution from Fe I lines remains large, while that derived from Fe II lines 
is narrow and compatible with no iron spread. We demonstrate that the metal-poor component claimed by Sim¬ 
merer et al. is composed by asymptotic giant branch stars that could be affected by non local thermodynamical 
equilibrium effects driven by iron overionization. This leads to a decrease of the Fe I abundance, while leaving 
the Fe II abundance unaltered. A similar finding has been already found in asymptotic giant branch stars of the 
globular clusters M5 and 47 Tucanae. We conclude that NGC 3201 is a normal cluster, with no evidence of 
intrinsic iron spread. 

Subject headings: stars: abundances — techniques: spectroscopic — globular clusters: individual (NGC 3201) 


1. INTRODUCTION 

The striking homogeneity in terms of iron abundance is 
usually considered as the main chemical signature of glob¬ 
ular clusters (GCs), indicating that these stellar systems were 
not able to reta in the gas ejected by supernova (SN) explo¬ 
sions (see e.g. [Renzini & Buzzoni 1986; Willma n & Straderl 
1201 2l) . For four decades, only one GC-like system, namely 
uj Centauri, was known to display an intrinsic dispersion 
in Fe ( see e.g.lF reeman & Rodge rslll975l:lNorris & Da Costal 
19951 lOriglia et alJ 120031 : iJohnson & Pilachowskil 120101 : 
Pancino et alJ 120 iTF) and this evidence brought to classify 
the system as the remnant core of a tidally disrupted dwarf 
galaxy accreted by the Milky Way. In the last few years, 
deep and extensive spectroscopic and photometric investiga¬ 
tions have revealed a more complex picture, with the discov¬ 
ery of other (massive) GCs harboring distinct sub-populations 
with different iron abundance , as Terzan 5 dFerraro et all 
l2009t lOrigiia et alJl201 ll l20ll iMassari et alJl2014l) and M2 
(lYong et al.ll2014jk or with large but uni-modal iron distribu¬ 
tions, as M54 dCarretta et akll2010al) and M22 (IMarino et alJ 
120091) . Further suggestions of intrinsic iron spread are only 
tentative at the moment, either because based on Ca II 
lines (NGC5824. lDa Costa. Held & Savianel2014l) or because 
of still disputed res ults (NGC 1851, ICarretta et alJ 1201 Obt 
IVillanova et al.ll2010l) . 

Among these candidate anomalous GCs, the case of 
NGC 3201 is controversial because different analyses pro¬ 
vide conflicting result^ about its level of iron homogene¬ 
ity. iGonzalez & Wallersteinl (i 19981 ) first analysed a sample 
of CTIO high-resolution spectra of 18 cluster stars, finding 
large iron variations (A [Fe/H] — 0.4 dex). However, an evi¬ 
dent trend between [Fe/H] and the effective temperature casts 
doubts about the reliability of their [Fe/H] distribution. The 
re-analysis of six of their targets performed by ICovev et~akl 
(120031) does not provide addition al clues. 

Further analyses by ICarretta et al.1 ( 2009 ) and 
iMunoz, Geisler & Villanoval (12013|) . based on 

high-resolution, high signal-to-noise ratio spectra 
(FLAMES @VLT and MIKE@Magellan, respectively), 


do not highlight similar spreads, ruling out large star-to- 
star variations. On the other hand, iSimmerer et all (120131) 
analysed UVES@FLAMES and MIKE@Magellan high- 
resolution spectra of 24 giant stars of the cluster, revealing a 
metallicity distribution as large as 0.4 dex (not explainable 
within the uncertainties) and with an evident metal-poor 
tail (5 out 24 stars). This iron_spread, qua litati vely sim¬ 
ilar to that observed in M22 (iMarino et alJ 120091 ). would 
make NGC 3201 the least massiv e GC (~ 1.1 x 10 5 M Q ; 
iMcLaughlin & van der MareH 12005 ) with evidence of SN 
ejecta retention. 

Recently, lLapenna et al] (120141) discovered effects due 
to the departure from local thermodynamical equilibrium 
(NLTE) in a sample of asymptotic giant branch (AGB) stars 
in the globular cluster 47 Tucanae. Such NLTE effects affect 
the abundances derived from Fe I lines (bringing to an under¬ 
estimate of 0.1-0.2 dex in [Fe/H]), but leave the abundances 
from Fe II lines unaltered. On the other hand, red giant branch 
(RGB) stars do not exhibit similar effects and the abundances 
from Fe I and Fe II nicely match each other. 

Based on this finding, it is well reasonable to ask 
whether the star-to-star scatter measured in the Fe content of 
NGC 3201 is genuine or it is due to such a spurious effect, be¬ 
cause of the inclusion in the sample of some AGB stars. Note 
that the NLTE effects can be easily unveiled by comparing the 
abundances derived from Fe I and Fe II separately. 

In this paper we re-analyse the spectra of the stars dis¬ 
cussed i n [Simmerer et al.l (120131) in light of the finding by 
lLapenna et alJ (120141) . We limit the analysis to the UVES- 
FLAMES sample, including also the 5 stars populating the 
metal-poor tail of the ISimmerer et all (120131 ) distribution. 

2. OBSERVATIONS 

High-resoluti on spec tra taken with UVES-FLAMES @VLT 
(iPasauini et alJ l2000t) for 21 giant stars members of 
NGC 3201 have been retrieved from the ESO archive. The 
spectra have been acquired with the UVES grating 580 
Red Arm CD#3, that provides a high spectral resolution 
(R—45000) and a large spectral coverage (—4800-6800 A). 
The spectra have been reduced using the dedicated ESO 
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pipeline^], performing bias subtraction, flat-fielding, wave¬ 
length calibration, spectral extraction and order merging. In 
each exposure one fiber is dedicated to sample the sky back¬ 
ground and used to subtract this contribution from each indi¬ 
vidual spectrum. 

Spectroscopic targets have been identified in our photomet¬ 
ric catalog, obtained by combining high resolution images 
acquired with the HST-ACS camera and wide-field images 
acquired with the ESO-WFI imager. Both the photometric 
datasets have been obtained through the V and I filters. A to¬ 
tal of 13 targets lie in the innermost cluster region, covered 
by ACS, while 8 stars are in the external region, covered by 
WFI. The membership of all the targets is confirmed by their 
very high radial velocity (< RVheiio >= +494.6 ± 0.8 km s _1 , 
a = 3.6 km s -1 ) that allows to easily distinguish the cluster 
members from the surrounding field stars. 

The position of the targets in the color-magnitude dia¬ 
grams (CMDs) is shown in Fig. Q] These CMDs have been 
corrected for differential reddening using the method de¬ 
scribed in Massari et ak| (f2012j) and adopting the extinction 
law by ICardelli. Clavton & Mathis (119891) . In order to cal¬ 
culate guess values for the atmospheric parameters of the 
target stars, we fitted the CMDs with an appropriate theo¬ 
retical isochrone from the BaSTI dataset (Pietrmfemietak 
2006), computed with an age of 11 Gyr (Marin-Franch et al. 
2009), Z= 0.001 and a-enhanced chemical mixture, finding a 
color excess E(B-V)= 0.31 mag and a true distance modulus 
(m-M)o=13.35 mag. 

Table 1 lists the main information about the targets, 
by adopting the sam e identification numbers used by 
Simmerer et aD (120131) who adopted the original names by 
Coteet al.l(l1994i)~ 

3. ANALYSIS 

Iron abundances have been derived with the package 
GALAFI dMucciarelli et al.ll2013al) by matching the measured 
and theoretical equivalent widths (EWs). Model atmospheres 
have been calculated with the code ATLASS0 We selected 
Fe I and Fe II lines predicted to be unblended at the UVES res¬ 
olution and at the typical atmospheric parameters and metal- 
licity of the observed stars, through the careful inspection 
of synthet ic spe ctra calculated with the SYNTHE package 
(lSbordoneii2005l) . Atomic data of the transitions of interest 
are from the last release of the Kurucz/Castelli linelistQ. The 
final iron abundances are based on ~ 130-150 Fe I and ~ 15- 
20 Fe II fines. EWs have been measured with DAOSPEC 
dStetson & Pancin o| [2008[). run iteratively by means of the 
package 4DAd’l (lMucci are iTil 1201 3l) . EW, oscillator strength 
and excitation potential for all the measured transitions are 
fisted in Table 2 (available in its entirety in the online ver¬ 
sion). 

3.1. Analysis with spectroscopic gravities 

First, we performed a fully spectroscopic analysis, as done 
by lSimmerer et all d2013i) . in order to verify whether we ob¬ 
tain the same evidence of a metallicity dispersion. In this 
analysis the atmospheric parameters have been constrained as 
follows: (a) for the effective temperatures (T e s) we requested 

1 http://www.eso.org/sci/software/pipelines/ 

2 http://www.cosmic-lab.eu/gala/gala.php 

3 http ://wwwuser. oats .inaf. it/castelli/sources/atlas9codes .html 

4 http://wwwuser.oats.inaf.it/castelli/linelists.html 

5 http://www.cosmic-lab.eu/4dao/4dao.php 


that no trend exists between abundances and excitation poten¬ 
tial, (b) for the surface gravities (log g) we imposed that the 
same abundance is obtained (within the uncertainties) from 
Fe I and Fe II fines, (c) for the microturbulent velocity (v tur b) 
we requested that no trend exists between abundances from 
Fe I fines and the reduced fine strength. The derived values of 
Vturb are based on ~ 130-150 Fe I fines distributed over a large 
interval of reduced EWs, with log(EW/A) ranging between - 
5.6 and-4.7. 

We derived an average [Fe I/H]=-1.46±0.02 dex (cr= 0.10 
dex), with a distribution ranging from -1.62 dex to -1.27 dex. 
The [Fe I/H] and [Fe II/H] abundance distributions are shown 
in Fig. [2] as generalized histograms. This result well matches 
that obtained by [Simmerer et all (12013l) that find an average 
abundance [Fe/H]=-1.48±0.02 dex (cr= 0.11 dexf] with a 
comparable iron range (A [Fe/H] ~0.4 dex). This spectro- 
scopic analysis fully confirms the claim by iSimmerer et all 
(120131) : when analysed with atmospheric parameters derived 
following the constraints fisted above, the stars of NGC 3201 
reveal a clear star-to-star scatte r in the iron content. The 5 
stars labelled as metal-poor by ISimmerer et al.l (120131) . with 
[Fe/H] <-1.58 dex, are the most metal-poor also in our metal- 
licity distribution. 

3.2. Analysis with photometric gravities 

As pointed out bv lLapenna et al. (2014?) , possible NLTE ef¬ 
fects in AGB stars can be easily detected by assuming pho¬ 
tometric values for log g and measuring Fe I and Fe II in¬ 
dependently. In order to see whether this effect is present 
in the NGC 3201 data, we adopted the following procedure. 
T e ff has been derived spectroscopically, by imposing the ex¬ 
citation equilibrium as described above. Thanks to the high 
quality of the spectra (with S/N ratio per pixel higher than 
100) and the large number of Fe I fines distributed over a 
large range of excitation potentials, very accurate spectro¬ 
scopic Teg can be estimated, with internal uncertainties of 
about 20-30 K. As a guess value we adopted 7A- calculated 
from the (V - I)q - T e ff calibration bv I Alonso et all (119991) and 
assuming a color excess E(B-V)=0.31 mag. Gravities have 
been derived from the Stefan-Boltzmann equation, assuming 
E(B-V)=0.31 mag, Q«- M)q= 1 3.35 mag, bolometric correc¬ 
tions from I Alonso et alJ d 19991) and a mass of 0.82 M 0 (ac¬ 
cording to the best-fit isochrone, the latter is a suitable value 
for RGB stars brighter than the RGB Bump magnitude level). 
Because T e ff is derived spectroscopically, the gravity is recom¬ 
puted through the Stefan-Boltzmann equation in each itera¬ 
tion according to the new value of T e ff. The mass value of 
0.82 Mq is appropriate for RGB stars but probably too high 
for AGB sta rs, bec ause of mass loss pheno mena during the 
RGB phase (iRoodll 1 973b fOri trlia et al]l2(Ml2(Ml20T4l) . In 

6 Note that iSimmerer et all <20131) adopted as solar reference value 
7.56 obtained from thei r own solar analysis, while we used 7.50 by 
[Grevesse_& Sauva l (19981) . Througout the paper we refer to the abundances 
bv lSimmerer*era!Tl2013l) corrected for the different solar zero-point. 

7 Recently, IJohnson et all 12014) analysed a sample of 35 AGB stars in 
47 Tucanae, finding no clear evidence of NLTE effects. Even if a detailed 
comparison between the two analyses is not the scope of this paper, we high¬ 
light som e main differences be tween the two works: the spectral resolution 
(48000 in lLapenna et all (20141) and 22 000 in IJohnson etalJ jOljlLjhe num¬ 
ber of Fe II lines (13 in lLapenna et all (2QI4| } and 4 in ljohnson e t al JI20 1 41 . 
on ave rage) and the adopted linelists for AGB and RGB stars I Laoe nna et all 
(2014D adopted the s ame linelist for both the groups of stars, at variance with 
lohnson et al. 1 2014 ) that used a linelist c onsistent, but not exactly the same, 
with that adopted by Cordero et al. (2014) where the reference RGB stars are 
discussed). 
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Fig. 1.— CMDs for central and external regions of NGC 3201 ( left and ri g ht panels, re spectively), corrected for differential reddening. Large circles are the 
spectroscopic targets flagged as metal-rich ([Fe/H]>-1.58 dex) in ISimmerer et alJ (120131) . while large triangles those identified as metal-poor ([Fe/H]<—1.58 
dex). 


TABLE 1 


Star 

RA 

Dec 

V 

I 

RVhelio 

T e ff 

logg 

Vturb 

[Fe I/H] 

[Fe II/H] 


(J2000) 

(J2000) 



(km/s) 

(K) 


(km/s) 

(dex) 

(dex) 


o3 

154.3084680 

-46.4125790 

13.76 

12.42 

495.9±0.6 

4730 

1.54 

1.45 

-1.27=1=0.03 

-1.35i0.05 

89 

154.3346190 

-46.3860090 

13.87 

12.56 

498.8T0.5 

4855 

1.63 

1.60 

-1.57=1=0.02 

-1.35i0.03 

91 

154.3359600 

-46.3355290 

13.81 

12.41 

490.44=0.6 

4705 

1.56 

1.50 

-I.36i0.05 

-1.35i0.05 

105 

154.3439854 

-46.4201000 

13.94 

12.59 

493.94=0.9 

4760 

1.61 

1.45 

-1.31i0.03 

-1.29i0.05 

124 

154.3601786 

-46.4140432 

12.85 

11.32 

498.54=0.2 

4375 

0.96 

1.55 

-1.37±0.02 

-1.40i0.05 

129 

154.3623900 

-46.4312480 

13.54 

12.13 

491.84=0.4 

4580 

1.36 

1.50 

-1.48i0.02 

-1.43i0.04 

181 

154.3876490 

-46.4126379 

13.89 

12.58 

497.54=0.4 

4920 

1.65 

1.65 

-1.62i0.03 

-1.37i0.04 

200 

154.3949420 

-46.3972062 

12.79 

11.30 

495.14=0.4 

4515 

0.99 

1.85 

-1.50i0.02 

-1.40i0.04 

222 

154.4031199 

-46.4257714 

12.69 

11.15 

491.04=0.5 

4355 

0.88 

1.65 

-1.47i0.02 

-1.45i0.05 

231 

154.4067690 

-46.4012923 

13.52 

12.18 

491.24=0.5 

4785 

1.45 

1.60 

-1.54i0.02 

-1.39i0.03 

240 

154.4091582 

-46.4277707 

13.91 

12.59 

495.24=0.6 

4855 

1.64 

1.55 

-1.61±0.03 

-1.38i0.04 

244 

154.4097946 

-46.4024205 

13.77 

12.38 

492.24=0.3 

4690 

1.51 

1.50 

-1.44i0.02 

-1.44i0.05 

249 

154.4104867 

-46.4270988 

12.99 

11.54 

498.14=0.3 

4545 

1.11 

1.60 

-1.56i0.02 

-1.42i0.04 

277 

154.4198397 

-46.4124756 

13.60 

12.18 

497.64=0.3 

4620 

1.39 

1.50 

-1.40i0.02 

-1.42i0.05 

279 

154.4208433 

-46.3954964 

13.34 

11.89 

499.04=0.3 

4555 

1.26 

1.50 

-1.43i0.03 

-1.43i0.04 

303 

154.4288640 

-46.3800960 

13.58 

12.20 

491.74=0.4 

4630 

1.39 

1.50 

-1.45i0.02 

-1.45i0.04 

308 

154.4304349 

-46.4108231 

13.76 

12.37 

487.24=0.4 

4620 

1.45 

1.45 

-1.50i0.02 

-1.45i0.04 

312 

154.4342122 

-46.4240651 

12.71 

11.10 

494.94=0.3 

4330 

0.85 

1.65 

-1.44i0.03 

-1.44i0.06 

332 

154.4478331 

-46.3982898 

12.96 

11.47 

499.74=0.5 

4500 

1.08 

1.65 

-1.54i0.02 

-1.38i0.04 

344 

154.4544130 

-46.4162830 

13.82 

12.41 

490.74=0.7 

4655 

1.49 

1.40 

-1.31i0.03 

-1.29i0.05 

374 

154.5004800 

-46.5194470 

13.64 

12.09 

498.74=0.6 

4595 

1.38 

1.50 

-1.41i0.03 

-1.36i0.05 










<[Fe I/H] > 

<[Fe II/H] > 










-1.46i0.02 

~1.40i0.01 


NOTE. — Main information of the target stars. Identification numbers are the same adopted bv ISimmerer et alJ d2013D . [Fe 
I/H] and [Fe II/H] have been obtained adopting photometric gravities. 


fact, iGratton et alJ 1201 Ob provide the masses for a sample of 
HB stars in NGC 3201, finding values between 0.62 and 0.71 
Mq. We initially analysed all the targets assuming the mass 
of a RGB star. Then, the AGB candidates, selected according 
to their position in 7 d1 -log g plane (as discussed in Section 
4) have been re-analysed by assuming the median value (0.68 
Mq) of the HB stars estimated by IGratton et all (12010(> . 

This method allows us to take advantage of the high-quality 
of the spectra, deriving accurate 7/n thanks to the large num¬ 
ber of transitions spanning a large range of excitation poten¬ 
tials. On the other hand, this approach does not require a 
fully spectroscopic determination of log g that is instead cal¬ 
culated using both photometric information and spectroscopic 


7^ff, avoiding any possible bias related to NLTE effects. 

The [Fe I/H] and [Fe II/H] abundances obtained with this 
method are listed in Table 1. The right panel of Fig.[2]shows 
the [Fe I/H] and [Fe II/H] distributions represented as gener¬ 
alized histograms obtained from this analysis. The two distri¬ 
butions turn out to be quite different. The iron distribution ob¬ 
tained from Fe I lines resembles that obtained with the spec¬ 
troscopic parameters (left panel of Fig. 0, with an average 
value of [Fe I/H]=-1.46±0.02 dex (o= 0.10 dex), while the 
distribution obtained from Fe II lines has a narrow gaussian- 
shape ([Fe II/H]>=-1.40±0.01 dex, o= 0.05 dex) pointing to 
a quite homogeneous iron content. 


























4 




[Fe/H] [Fe/H] 


FlO. 2.— Generalized histograms for [Fe I/H] (empty histogram) and [Fe 
gravities (left panel) and with photometric gravities (right panel). 


TABLE 2 

Star identification number, wavelength, 

OSCILLATOR STRENGTH, EXCITATION POTENTIAL 
AND MEASURED EWS FOR ALL THE USED 
TRANSITIONS. 


Star A Ion log(gf) E.P. EW 
(A) (eV) (mA) 


63 

4791.246 

Fel 

-2.435 

3.270 

26140 

63 

4834.507 

Fel 

-3.330 

2.420 

27.90 

63 

4842.788 

Fel 

-1.530 

4.100 

16.30 

63 

4892.859 

Fel 

-1.290 

4.220 

29.50 

63 

4911.779 

Fel 

-1.760 

3.930 

19.50 

63 

4917.230 

Fel 

-1.160 

4.190 

35.90 

63 

4918.013 

Fel 

-1.340 

4.230 

28.30 

63 

4950.106 

Fel 

-1.670 

3.420 

54.10 

63 

4962.572 

Fel 

-1.182 

4.180 

32.30 

63 

4969.917 

Fel 

-0.710 

4.220 

46.00 

63 

4985.253 

Fel 

-0.560 

3.930 

75.30 

63 

5002.793 

Fel 

-1.530 

3.400 

67.20 

63 

5014.942 

Fel 

-0.303 

3.940 

86.70 

63 

5022.236 

Fel 

-0.560 

3.980 

76.40 

63 

5028.126 

Fel 

-1.123 

3.570 

68.30 


NOTE. — This table is available in its entirety 
in a machine-readable form in the online journal. A 
portion is shown here for guidance regarding its form 
and content. 

3.3. Uncertcinties 

Internal uncertainties in the derived Fe abundances have 
been calculated by adding in quadrature two sources of un¬ 
certainties: 

(1) those arising from the EW measurement. For each target, 
we estimated this term as the line-to-line dispersion normal¬ 
ized to the root mean square of the number of lines. Because 
of the high quality of the used spectra, the line-to-line scatters 
are smaller than 0.1 dex, leading to internal uncertainties of 
about 0.005-0.008 dex for [Fe I/H] and of about 0.010-0.025 
dex for [Fe II/H]. 

(2) those arising from the atmospheric parameters. To 
estimate this term, we follow the approach described by 
ICavrel et al.l (120041) to take into account the covariance terms 
due to the correlations among the atmospheric parame- 


(grey histogram) obtained from the analysis performed with spectroscopic 


TABLE 3 

Abundance uncertanties due to the 

ATMOSPHERIC PARAMETERS FOR THE STARS #63 AND 
#89. 


Ion 

Parameters 

Uncertainty 

(dex) 

<5 Tcff 
±50 K 
(dex) 

<51ogg 

±0.1 

(dex) 

^ V turb 

dzO.l km/s 
(dex) 


#63 (RGB) 

Fe I 

±0.03 

±0.05 

±0.00 

4=0.03 

Fell 

±0.04 

4=0.02 

±0.04 

±0.03 

#89 (AGB) 

Fe I 

±0.02 

±0.04 

±0.00 

4=0.03 

Fell 

±0.03 

4=0.02 

±0.04 

±0.02 


NOTE. — The second column is the total un¬ 

certainty calculated according to ICavrel et all 12004 ). The 
other columns list the abundance variations related to the 
variation of only one parameter. 

ters. For each target, the temperature has been varied by 
± 1 CT 7 df , the gravity has been re-calculated through the Stefan- 
Boltzmann equation adopting the new values of T c a and the 
microturbulent velocity re-optmized spectroscopically. 

Table 1 lists the total uncertainty including both the terms 
(1) and (2). Also, Table 3 shows for two representative tar¬ 
gets (one RGB and one AGB star) the abundance uncertainty 
obtained following the prescriptions by ICavrel et al.l (120041) 
(second column) and those obtained with the usual method of 
independently varying each parameter (an approach that ob¬ 
viously does not take into account the correlation among the 
parameters and can over-estimate the total uncertainty). 

4. DISCUSSION 

In this paper we present a new analysis of the UVES- 
FLAMES spectra of 21 member stars of NGC 3201 already 
discussed in ISimmerer et al.l (120131) . The Fe abundances 
have been calculated both using spectroscopic gravities (ob¬ 
tained by imposing the ionization balance between Fe I and 
Fe II abundances) and photometric ones (obtained through the 
Stefan-Boltzmann equation). The two methods provide dif¬ 
ferent results concerning [Fe I/H] and [Fe II/H]. In particular, 
the use of spectroscopic log g provides a wide [Fe/H] distri¬ 
bution (al large as ~=0.4 dex), in agreement with the finding of 


II/H] 
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ISimmerer et al.l(l2013l) who adopted the same method. On the 
other hand, when photometric gravities are used, the [Fe I/H] 
distribution remains quite large, while that of [Fe II/H] is nar¬ 
row. We compute the intrinsic spread of the two Fe distribu¬ 
ti ons adopting the Maximu m Likelihood algorithm described 
in iMucciarelli et al.l ( 2012 ). Concerning [Fe I/H] we derive 
ai n , = 0.09 ±0.01 dex, while for [Fe II/H] o int = 0.00 ±0.02 
dex. Hence, the [Fe II/H] distribution is compatible with no 
iron spread. _ 

ISimmerer et al. : feol3i) highlight that the 5 most metal-poor 
stars of their sample are bluer than the other stars, as expected 
in cases of a lower metallicity. The left panel of Fig.[3]shows 
the position of the targets in the T e ff-log g plane, with su¬ 
perimposed, as reference, two isochrones with the same age 
but different metallicity: Z= 0.001 (solid line) and Z= 0.0006 
(dashed line). The RGB of the isochrone with Z= 0.0006 over¬ 
laps the position of the AGB of the Z= 0.001 isochrone. Seven 
targets (including the 5 candidate metal-poor stars) are located 
in a position compatible with both the scenarios: metal-poorer 
RGB or AGB at the cluster metallicity. These seven targets 
have average abundances of [Fe I/H]=-1.57±0.01 dex and [Fe 
II/H]=-1.41±0.01 dex, while the RGB stars have [Fe I/H]=- 
1.42±0.02 dex and [Fe II/H]=-1.40±0.01 dex. 

However, if these stars were metal-poor RGB stars, the iron 
abundance derived from Fe I lines should be in agreement 
with that obtained from Fe II, since NLTE effects are not ob¬ 
served in RGB stars of comparable luminosity and metallicity. 
The right panel of Fig. [3] shows the behaviour of ([Fe I/H] - 
[Fe II/H]) as a function of the difference between the spectro¬ 
scopic values of 7±f and those obtained from the projection 
along the RGB of the best-fit isochrone. The stars located 
along the RGB (A7/if ^0) have similar [Fe I/H] and [Fe II/H] 
abundances, compatible with no NLTE effects. The stars hot¬ 
ter than the reference RGB have differences between [Fe I/H] 
and [Fe II/H] ranging from -0.07 to -0.22 dex, with a mean 
value of -0.15 dex. The grey region marks the expected po¬ 
sition for metal-poorer RGB stars: they should be hotter than 
the reference RGB, but with [Fe I/H] - [Fe II/H]~0, as com¬ 
monly measured in the RGB stars. This reveals the true na¬ 
ture of these stars: they are genuine AGB stars, with the same 
metallicity of the cluster (as measured from their Fe II lines) 
but affected by NLTE effects leading to a systematic decrease 
of [Fe I/H] . This is the same e ffect observed by Iva ns et all 
(2001) and lLaoenna et al.l (120141) in the AGB stars of M5 and 
47 Tucanae, respectively. 

A direct inspection of the spectra reveals the different be¬ 
haviour of Fe I and Fe II lines in AGB and RGB stars. Fig. |4] 
shows three Fe I lines (chosen with different excitation po¬ 
tential) and one Fe II line in the spectra of the AGB star 
#89 (upper panels) and of the RGB star #303 (lower panels). 
Synthetic spectra calculated with the appropriate atmospheric 
parameters and the metallicity derived from Fe II lines (red 
lines). In the upper panels we also show the synthetic spec¬ 
trum computed with the average [Fe I/H] (blue dashed line). 
Clearly, the synthetic spectrum assuming the [Fe II/H] abun¬ 
dance well reproduces all the observed lines in the case of 
the RGB star, while it fails to fit the Fe I lines observed in 
the AGB star, regardless of the excitation potential (pointing 
out that this effect cannot be attributed to inadequacies in the 
adopted T e g). On the other hand, the abundance derived from 
Fe I lines is too low to well reproduce the depth of the AGB 
Fe II line plotted in Fig. [4] This clearly demonstrates a differ¬ 
ent behaviour of iron lines in AGB and RGB stars. 


It is worth noting that this behaviour is somewhat puzzling, 
because theoretical models do not predict significant differ¬ 
ences in the NLTE corrections for stars in the parameter space 
covered by our targe ts. For instance, the grid of NLTE cor¬ 
rection s computed by Bergem ann et alJ (120121) and lLind et alJ 
(120121) predicts that the Fe I lines in AGB and RGB stars 
should be affected in a very similar way at the metallicity of 
NGC 3201. However, some additional effects/mechanisms 
could play a role in the AGB photospheres, leading to the de¬ 
parture from the LTE condition, which are not yet accounted 
for in the available theoretical calculations. 

We checked whether the mass assumed for AGB stars could 
change our conclusions. The mass distribution of HB stars in 
NGC 3201 provided bv lGratton et al.l (120101) ranges from 0.62 
to 0.71 M 0 , with a median value of 0.68 Mq. The minimum 
and maximum mass values correspond to a difference in log g 
of 0.06, leading to a variation in [Fe II/H] of only 0.02-0.025 
dex (and a variation in [Fe I/H] of 0.002-0.005 dex). On the 
other hand, if we adopt the RGB mass (0.82 Mq ) for all tar¬ 
gets, as it is often done because of the difficulty to observa- 
tionally distinguish between RGB and AGB stars, the values 
of [Fe II/H] increase by only 0.03 dex with respect to esti¬ 
mate obtained assuming 0.68 Mq. Hence, the precise value 
of the adopted mass (within a reasonable mass range) cannot 
reconcile the difference between [Fe I /H] and [Fe II/H]. 

As additional check, for each AGB stars, the stellar mass 
has been varied until the ionization equilibrium was satis¬ 
fied. The derived values range from ~0.2 and ~0.5 M 0 : such 
masses are too low with respect to the mass distribution of the 
HB stars derived bv lGratton et al.l (120 10i) . In particular, for the 
stars #89, #181 and #240, that exhibit the largest difference 
between [Fe I/H] and [Fe II/H] (~-0.2 dex), a satisfying ion¬ 
ization equilibrium can be reached only with masses smaller 
than 0.2-0.25 Mq, which are very unlikely values for globular 
cluster AGB stars. 

5. CONCLUSIONS 

We demonstrated that the observed intrinsic star-to-star Fe 
scatter in the GC NGC 3201 is due to unaccounted NLTE 
effects in the spectroscopic analysis of some AGB stars in¬ 
cluded in the sample. These stars suffer from NLTE effects 
driven by the iron overionization, a mechanism that affects 
mainly the less abundant species like Fe I, but has no a sig¬ 
nificant effect on the dominant species (e.g. Fe II). When the 
gravity of these stars is obtained spectroscopically, forcing to 
have the same abundance from Fe I and Fe II lines, the derived 
[Fe/H] abundance turns out to be under-estimated. 

Our findings confirm the conclusion by iLanenna et alJ 
(120141) that the chemical analysis of samples of stars including 
both AGB and RGB stars, and based on spectroscopic gravi¬ 
ties, can lead to spurious broadening of the iron distribution. 

We conclude that NGC 3201 is a normal GC, without evi¬ 
dence of intrinsic iron scatter. In light of this result, it is not 
necessary to suppose that NGC 3201 was more mass ive in 
the pa st to retain the SN ejecta, as invoked bv ISimmerer et alJ 
(120131) . 

We warmly thank the anonymous referee for suggestions 
that helped improving the paper. This research is part of the 
project COSMIC-LAB (http://www.cosmic-lab.eu) funded by 
the European Research Council (under contract ERC-2010- 
AdG-267675). 
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Fig. 3.— Left panel : position of the targets in the T^-log g plane. A BaSTI isochrone with age of 11 Gyr, Z=0.001 and a-enhancement chemical mixture 
(solid line) is superimposed for sake of comparison. The dashed line indicates the position of the RGB for a BaSTI isochrone with age of 11 Gyr an Z=0.0006. 
Grey triangles are the metal-poor stars of ISimmerer et al.l | |2013 ). Empty grey squares are the stars in our analysis with [Fe I/H] - [Fe II/H] <—0.1 dex. Right 
panel: Behaviour of [Fe I/H] - [Fe II/H] as a function of the difference between the spectroscopic temperature of each target and the photometric value estimate 
from the isochrone with Z=0.001 shown in the left panel. Same symbols as in the left panel. The grey light area indicates the expected mean locus for metal-poor 
([Fe/H]~-1.6 dex) RGB stars. 
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Fig. 4.— Spectral regions around three Fe I lines with different excitation potential and one Fe II line, for the AGB star #89 (upper panels) and the RGB star 
#303 (lower panels). Synthetic spectra calculated with the corresponding atmospheric parameters (see Table 1) and adopting the average iron abundance derived 
from Fe II lines are superimposed as red curves. The blue dashed curve shown in the upper panels is the synthetic spectrum calculated with the iron abundance 
derived from Fe I lines. 











